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Abstract 

The phase diagram of Yukawa particles confined between two parahel hard wahs is calculated at 
zero-temperature beyond the bilayer regime by lattice-sum-minimization. Tuning the screening, a 
rich phase behavior is found in the regime bounded by stable two-triangular layers and 3-square 
layers. In this regime, alternating prism phases with square and triangular basis, structures derived 
from a hep bulk lattice, and a structure with two outer layers and two inner staggered rectangular 
layers, reminiscent of a Belgian waffle iron, are stable. These structures are verifiable in experiments 
on charged colloidal suspensions and dusty plasma sheets. 

PACS numbers: 82.70.Dd, 64.70.K- 
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If a system is confined to a thin slit geometry, its properties are drastically altered with 



respect to the bulk 



,|2j. In particular, the freezing transition in confined geometry occurs 
into multilayered crystals whose structure is dictated by an interplay of the interparticle 
interaction and the confining potential . For hard spheres between two parallel hard 

plates of finite separation D, different solid lattices are getting stable with the following 
cascade for increasing width D: li D equals the hard sphere diameter, a crystalline monolayer 
is stable which possesses a two-dimensional triangular (A) lattice symmetry. This layer 
buckles upon increasing h such that a two-layer situation is realized. For higher 

h, there are stable structures that correspond to two intersecting square lattices (2n) or 

nn fin 

triangular lattices (2 A) [9|, Il0| with a two-layer rhombic phase in between [6|, Il0|. Beyond 
bilayers, the transition from the 2A to the 30 (simple quadratic trila yer ) phase are mediated 
by prism phases which are made up of alternating prism-like arrays ll|, II4]. In experiments 



on strongly screened colloidal suspensions, however, crystals derived from the hexagonal- 
closed-packed geometry were found as interpolating structures which correspond to four 
or more layers jisl . 14\. Beyond the SD phase, the structures are even more complicated 



]J,ll5|,ll6|. 



Another situation concerns the confined Wigner crystal of repulsive point charges. This 
classical Coulomb system has been studied when confined between two plates and a bilayer 



structure is always stable 



17 1 even for high separation distances h [18|. Here we consider 



a Yukawa system interpolating between the Coulomb interaction and hard spheres 19|. In 
this case, there are two possibilities for the rhombic bilayer 201]. 
can be thought of as charged colloids [2l| or dusty plasma sheets 



Ixperimental realizations 



22j. 



In this paper we study the multilayering for a Yukawa system confined between two hard 



walls 



23| interpolating between the Coulomb and hard-sphere limit. For low screening, ap- 



propriate to dusty plasmas, and an external parabolic potential, multilayerin g w as examined 
by computer simulation and a simple sheet model by Totsuji and coworkers [24?]. Structural 
details of the multilayers, however, were neglected in this shell model approach. Here we 
focus on the range of plate distances between the 2A to SD regime and resolve the fine 
structure using up to four layers as a possible candidate. Apart from prism phases and 
structures derived from the hexagonal-closed-packed geometry, we find a stable four-layer 
crystal with a structure resembling a Belgian waffle iron {BWT). This 5VI/7-structure which 
consists of two outer layers and two inner staggered rectangular crystals does not exist for 
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FIG. 1: Zero temperature phase diagram of a Yukawa system confined between two hard plates with 
screening strength A and reduced density rj, in the multilayer regime between 2 triangular and 3 
square layers. Stability regions of the Belgian waffle iron (BWI) phase, prism phases with triangular 
(2-Pa) and square bases (2Pg), 2hcp-like and 2hcp(100) phases are shown. For comparison, the 
hard-sphere limit of infinite screening is shown separately in the top. The arrow indicates a path 
of constant screening with increasing density. 



hard spheres and unscreened charges but is stable in the intermediate screening regime. A 
control over the achievable crystalline multilayered structures has important applications to 



25|. 



the fabrication of nanosieves and filters with desired porosity 

In our model, we consider N classical point-like particles interacting via the Yukawa 
pair-potential 

V{r) = Vo^^, (1) 
Kr 

where r is the interparticle distance, 1/k the screening length, and Vq denotes an energy 
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amplitude. The particles are confined between two parallel hard plates of area A and sep- 
aration D which is taken along the z-direction for convenience. At zero temperature, for a 
given reduced area density rj = ND'^/A, the system will minimize its total potential energy 
and the resulting optimal structure will be independent on the single energy scale Vq but 
will depend on the reduced inverse screening length A = kD. By varying this parameter 
A, one can interpolate between the unscreened Coulomb limit (A — 0) and the hard-sphere 
limit (A —>■ oo) where the interaction is getting very harsh. 

For a given screening parameter A and reduced density t] we have performed lattice sum 
minimizations of a set of candidates. While the stability regimes of mono- and bilayers have 
been addressed in previous work 20(], we focus here on the multilayer regime beyond the 
stability of the staggered triangular bilayer 2A phase. This can be achieved by increasing 
the reduced density rj. There is another region at even higher D where the trilayer phase 
Sn is stable. Our major goal is to explore the corresponding distance range between the 2A 
and Sn phase and check for the stability of intervening crystalline multilayers. 

As possible candidates for our lattice sum minimization problem, we consider three di- 
mensional crystals with a two dimensional periodicity in x and y-direction whose primi- 
tive cell is a parallelogram containing n particles. This parallelogram is spanned by the 
two lattice vectors a = a(l,0,0) and b = a7(cos^, sin^, 0), where 7 is the aspect ratio 
(7 = |b|/|a| = b/a) and 6 is the angle between a and b. Furthermore the n particles are 
distributed, not necessarily evenly, on m layers in the xy-plane. Hereby we restrict ourselves 
to layered situations with an up-down inversion symmetry in the averaged occupancy re- 
flecting the up-down symmetry of the confining slit. Under this sole restriction, we consider 
possible candidates with n = 2, ...,8 and m = 1, ...,4 up to symmetric four-layer structures 
with a basis of up to 8 particles. For given rj and A, the total potential energy per particle 
is minimized with respect to the particle coordinates of the basis and the cell geometry (7 
and 6). The resulting stability phase diagram is shown in fig. [T] 2^. 

First, the boundary of the stability domain of the staggered triangular phase 2 A and the 
staggered quadratic 30 are presented, see fig. [TJ For convenience, the structure of these two 
phases are sketched by a top and side view in fig. [2l In between, we observe various stable 
four-layer structures. For low screening (A < 4), however, there is a direct transition from 
2A to 3n without any intermediate multilayered crystal consistent with earlier simulations 
and theoretical results 2J]. The 2 A SD phase boundary moves to higher densities as the 
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FIG. 2: (a) and (b): Top view and side view of the 2A phase, (c) and (d): Top view and side view 
of the SD phase. Note that the layer distance does not necessarily coincide with the length of the 
in-plane square. Particles in different layers are drawn in different colors (green, red, blue). 



screening A is lowered so that there are only bilayers for the unscreened Coulomb system 
(A 0), see fig.[I]|30|. 

We now discuss the more complicated intermediate structures of fig. [U At finite A > 4, 
the 2A transforms discontinuously into a tetra-layered structure (m = 4) with n = 8. This 
novel structure shown in fig. [3] is characterized by two inner staggered rectangular lattices 
(see the red and yellow spheres in fig. [3]) and two outer layers with rhombic stripes which 
are centered relative to the inner rectangular lattice. In principle this structure can be 
obtained by a continuous transformation of the 2A phase. In analogy to Belgian waffles 
which possess an internal rectangular structure we call this double staggered situation a 
"Belgian waffle iron" {BWI) structure, as the corresponding iron mimics the outermost 
structure qualitatively [311]- The rhombus in the outermost layers is almost symmetric, i.e. 
the corresponding anisotropy as characterized by the ratio of width wi (see fig. [3]) and side 
length of the rhombic is close to 1 (within 1 percent) and only slowly grows with increasing 
density rj. 

Next, prism phases with a triangular basis shape 2Pa and a square basis shape 2Pn 
each with n = 6 and m = 4 (see fig. Hj), are stable. In detail, these phases consist of 
alternating double layered prism-like arrays with a triangular or square basis structure (see 
fig. HI). Within the stability regime of prism phases the distances W2 and ws indicated in fig. 
m change about few percents such that with increasing rj one can notice a slight decrease of 
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FIG. 3: Different views of the Belgian waffle iron [BWI) structure: (a) Top view of an outer layer 
and the next inner layer. The latter corresponds to a rectangular lattice (red particles) with a 
rhombic decoration (green particles), (b) Side view of all four layers, (c) Top view of all four 
layers. The width wi of the rhombus is shown in (a). 
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FIG. 4: Different views of the prism phases with triangular and quadratic shaped basis: (a) Top 
view of an outer layer and the next inner layer with the appropriate width W2 and (b) perspective 
side view of 2P^. (c) Top view of an outer layer and the next inner layer with the appropriate 
width W3 and (d) side view of 2Pq. The alternating prism like arrays are indicated with white 
lines in (b) and (d). 
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FIG. 5: Top views (a), (c) and side view (b) of 2hcp-like as well as top views (d), (f) and side view 
(e) of 2hcp(100). For 6^0, one recovers 2hcp(100). The angle 6 results from the additional shift 
in x-direction of the rectangular bilayers and is shown in (b) . 

these distances. 

Finally there are tetralayered structures derivable from the hep lattice as discussed re- 



cently in ref. 



13 



14l | which we accordingly call 2hcp-like and 2hcp(100) phases (see fig. Ej). 



The (lOO)-plane of the hexagonal-closed-package, possessing a rectangular shaped basis as 
sketched in fig. [5l can be thought of as two double layered rectangular arrays which lie on 
the top of each other but shifted to each other in the y-direction. The rectangular bilayers of 
the 2hcp-like structure, on the other hand, experience a further shift in x-direction, resulting 
in an angle 6 (see fig. El). 

For increasing screening constant A, the topology of the phase diagram reveals the fol- 
lowing cascades interpolating between 2A — > SD: 

2A 30, 
2A ^ BWI 30, 
2A ^ BWI 2Pa 2Pa 30, 
2A ^ BWI 2Pa 2hcp(100) 2Pa 2Pn 30 



(2) 
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with a remarkable reentrance of the 2Pa and finally 

2A ^ BWI 2Pa 2hcp-like ^ 2hcp(100) 

2Pa ^ 2Pn ^ 3n. (3) 

The hard-sphere limit involves the cascade 

2A ^ 2hcp-like ^ 2hcp(100) 2hcp-like 2Pu 30 

which is qualitatively different from those at finite screening as the BWI and the 2Pa are 
missing. Therefore even at high A ~ 100 there are considerable deviations from the ultimate 
hard-sphere limit. 

All transitions involved here for finite A are first order except the 2hcp-like — > 2hcp(100) 
one which is continuous. This is illustrated in fig. [6] where the distance Ac? between the 
outermost and the next inner layer is shown along the arrows indicated in fig. [T|, i.e. for 
increasing density t] at fixed screening A = 20 and A = 100 leading to the full transition 
cascade (see ([3])). The geometric meaning of Ad is sketched in the inset of fig. [61 The 
reduced quantity Ad/D can serve as an order parameter for multilayering in the sense, 
that Ad/D = 0.5 for centered trilayers and (formally) Ad/D = for symmetric bilayers. 
Multilayers beyond m = 3 are characterized by an intermediate value of Ad/D. Moreover 
for an equilayer-spacing, Ad/D is l/(m — 1). For A = 100, clearly, the transition from 
2A to BWI is accompanied with a jump from zero to Ad/D = 0.0221 {Ad/D = 0.0933 for 
A = 20) implying a first-order transition although a second-order transition would have been 
not precluded by lattice symmetry. The subsequent BWI —>■ 2Pa transition is discontinuous 
(though the jump in Ad/D is tiny in fig. [6]) while the 2hcp-like — > 2hcp(100) transition is 
continuous. The latter transition occurs as the angle 6 approaches zero as indicated in 
fig. [5l For A = 20, there are only discontinuous transitions along the cascade ([2]). In the 
limit of hard spheres, volume fraction calculations among the considered candidates show 
that the transitions 2A 2hcp-like — > 2hcp(100) 2hcp-like are of second order, while 
the transition 2hcp-like 2Pn is of first order. The transition 2Pn 30 occurs in the 
hard-sphere limit continuously, as opposed to finite A. 

To summarize, we investigated the stability phase diagram of a confined classical one- 
component Yukawa system in the transition regime 2A 30. Various crystalline tetra- 
layer structures are intervening including a novel one with a " Belgian waffie iron structure" 
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FIG. 6: Reduced distance Ad/D between the outermost and the next inner layer as a function of 
density tj for fixed A = 20 and A = 100. The limits of a trilayer and bilayer situations are indicated 
by the dashed lines. The geometric meaning of Ad is shown in the inset, developing of the distance 
between inner and outer layers. 

(BWI), phases of alternating prisms 2P^, 2Pn, and structures derived from a bulk hep 
structure (2hcp-like and 2hcp(100)). For the unscreened Coulomb system only bilayers are 
stable while in the opposite (hard-sphere) limit the sequence for increasing density is: 2A — >■ 
2hcp-like 2hcp(100) 2hcp-like 2Pn — * 30 [13]. Our theoretical results are verifiable 
in experiments on colloidal suspensions. In fact, all phases were found in recent experimental 
investigations except the BWI structure which, however, resembles experimentally deter- 



mined superstructures 



21| . Furthermore strongly interacting dust particles in plasmas are 



typically in the intermediate screening regime at low reduced temperatures and are therefore 
other candidates for an experimental system where the multilayer phases proposed here are 
in principle observable. 

Future work should focus on the inclusion of higher layer numbers m. We have examined 
the special case with m = 6 layers which was called pre-3n in ref. 13|. This structure. 



however, was not found to be stable in the range of A considered in this work. Next, charged 



walls will lead to a soft exponentially-screened wall-particle interactions 



32j rather than a 



hard wall which would introduce more parameters in the corresponding model 24] • this 



situation, the multilayering scenario is expected to 



the inner layers. The effect of finite temperature 33| needs more investigation; temperature 



avor phases which have more weight in 



can be viewed as an additional coordinate in the two-dimensional phase diagram presented 
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in fig. [U We expect, however, no drastic chan ge i n the phase diagram topology. Finally we 



34| might be a versatile tool to explore even 



remark that genetic minimization algorithms 
more complicated crystalline layer structures. 
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